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^ ■ Abstract 

O 

m . 

0^ . The infinite-C/ Anderson model is applied to transport through a quantum 



dot. The current and density of states are obtained via the non-crossing ap- 



ed 

. proximation for two spin-degenerate levels weakly coupled to two leads. At 

O ■ low temperatures, the Kondo peak in the equilibrium density of states strongly 



enhances the linear-response conductance. Application of a finite voltage bias 
reduces the conductance and splits the peak in the density of states. The 
split peaks, one at each chemical potential, are suppressed in amplitude by 
a finite dissipative lifetime. We estimate this lifetime perturbatively as the 
time to transfer an electron from the higher chemical potential lead to the 
lower chemical potential one. At zero magnetic field, the clearest signatures 
of the Kondo effect in transport through a quantum dot are the broadening, 
shift, and enhancement of the linear-response conductance peaks at low tem- 
peratures, and a peak in the nonlinear differential conductance around zero 
bias. 
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I. INTRODUCTION 



The Kondo effect has been a focus of condensed matter research for many years. 
Its essence - the crossover from weak to strong couphng between a magnetic im- 
purity and a conduction-electron sea as temperature is lowered - has inspired both 
theory and experiment. While most aspects of the problem are now well understood, 
studies have traditionally been confined to equilibrium properties.0 For the case of a 
magnetic atom embedded in a bulk metal, achieving nonequilibrium may be deter- 
ring, but it is not in the case of "artificial atoms" .i In particular, we predict that a 
quantum dot weakly coupled to its leads is a Kondo system in which nonequilibrium 
can be routinely achieved.lii More generally, an impurity or defect level in a small 
structure where the applied bias is dropped over a mesoscopic length!'! will be a 
nonequilibrium Kondo system. 

Anderson's mo dell for a Kondo impurity - a site with discrete, interacting lev- 
els coupled to a band - has already been used successfully to describe experiments 
on quantum dots.i^ The discrete spectrum of a single dot has been observed by 
transport0^ii and capacitance0 spectroscopy, while the strong on-site Coulomb 
interaction is recognized0 as the origin of periodic conductance oscillations.liiliiJii 
However, it is only the high temperature regime that has been explored experimen- 
tally, while it is at low temperatures that the Kondo effect emerges. 

Since the Anderson Hamiltonian describes the quantum dot, at low temperatures 
the dot must behave as a Kondo impurity. In fact, Glazman & Raikh0 and Ng & 
Leei have argued that at zero-temperature equilibrium the Kondo resonance in the 
density of states of spin-degenerate levels will produce perfect transparency of a 
quantum dot symmetrically coupled to its leads. More precisely, for all chemical 
potentials between eo and eo + U, where eo is the bare-level energy and U is the 
interaction energy (Fig. 0), the dot will have the conductance of an open channel, 
2e^/h. This is to be contrasted with the situation at temperatures larger than F, 
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the elastic width of the levels, where the conductance consists of two resonances, 
at eo and at eo + U. Since the chemical potential of a quantum dot can effectively 
be swept by changing the voltage on a nearby gate, the Kondo effect will have a 
striking experimental signature in low-temperature transport through a quantum 
dot. 

Until now, however, only qualitative predictions have been made for exper- 
imental observation of the Kondo effect in transport through a quantum dot. 
Specifically, raising the temperature above the relevant Kondo temperature is pre- 
dicted to suppress the peak in the density of states, and, consequently, reduce the 
conduct ance.0lli As the Kondo temperature near the conductance peak at eo de- 
pends exponentially on the chemical potential, T^(/i) ~ exp[— 7r(yU — eo)/r], one 
expects that the Kondo effect will enhance the conductance mainly for cq < ^ 
eo + few r. Accordingly, Ng and Leelll predicted that the finite temperature con- 
ductance vs. gate voltage will consist of pairs of asymmetric peaks,0 separated by 
the Coulomb-interaction energy U. In this work we present, for the first time, a 
quantitative calculation of the lineshape of these conductance peaks, via the non- 
crossing approximation.^ We show that at experimentally accessible temperatures 
the Kondo effect will leave the conductance peaks symmetric. The Kondo effect will 
manifest itself, instead, in the broadening of the peaks, the enhancement of their am- 
phtude, and the shift in their positions towards each other (for each spin-degenerate 
pair), as the temperature decreases. 

Furthermore, as the leads coupled to a quantum dot are easily biased to nonequi- 
librium, new physical questions which were not relevant to magnetic impurities can 
also raised. In particular, what happens to the Kondo effect out of equilibrium?i0 
Since at equilibrium the Kondo peak in the density of states occurs at the chemical 
potential, the presence in nonequilibrium of two chemical potentials must have a 
dramatic effect. 

In this paper we try to answer the question of what new behavior is present in 
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the Anderson model out of equilibrium, and to make quantitative predictions for 
experiment. Generalizing the non-crossing approximation to nonequilibrium, using 
the Keldysh formalism, we find that a voltage bias between the left and right leads 
causes the Kondo peak in the density of states to split, leaving a peak in the density 
of states at the chemical potential of each lead (Fig. |^). The amplitudes of these 
split peaks are suppressed by a finite nonequilibrium lifetime, due to dissipative 
transitions in which electrons are transferred from the higher chemical potential 
lead to the lower chemical potential one. The narrowness of the Kondo peak in 
the density of states, and its splitting and suppression, lead to a rapid decrease of 
conductance with increasing bias. The resulting peak in the nonlinear differential 
conductance is likely to be the most accessible experimental signature of the Kondo 
effect in quantum dots. 

We begin this paper with a general formulation of nonequilibrium transport 
through an Anderson impurity in the limit of an infinite on-site interaction energy, U. 
(section |11 Ap . Short discussions of the mapping into a "slave-boson" Hamiltonian 
and of the Keldysh formalism are presented. The non-crossing approximation is 
then introduced (section [lllj| ) and the numerical methods outlined (section |11 C| ). 
Results of the non-crossing approximation are presented for both equilibrium and 
nonequilibrium transport (section [II D|) . The theoretical interpretation of the results 



is discussed (section [III A|) as well as the implications for experiment (section piB[) . 
An appendix is included to demonstrate the current conserving property of the 
non-crossing approximation. 

II. THE NONEQUILIBRIUM ANDERSON MODEL 
A. General Formulation 
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1. The Model 

We model the quantum dot and its leads by the Anderson HamiltonianS 

H = ^kaCt^Cka + e^C+C^ + 1^^^ + iVkaC^^C^ + h.c), (1) 

a;k<^L,R f cr cr'jtcr ij;keL,R 

where c^^{cka) creates (destroys) an electron with momentum k and spin a in one 
of the two leads, and c+(Co-) creates (destroys) a spin-a electron on the quantum 
dot. The spin quantum number a may also represents orbital degeneracies as in 
the magnetic impurity problem,! though, experimentally, these degeneracies are 
likely to be lifted by disorder in quantum dots. In the following we will focus on 
spin-degenerate states. The third term describes the Coulomb interaction among 
electrons on the dot. We assume that f/ ^ oo, forbidding double occupancy. This 
is appropriate for quantum dots where, typically, f/(~ ImeV) is a hundred times 
larger than the coupling to the leads,!! r(~ lOyueV). The fourth term describes the 
hopping between the leads and the dot, and determines this coupling strength via 

r^(«'(^) = 2vr \Vk.?K^-ek.). (2) 

fceL(fl) 

Our aim is to calculate the current through the quantum dot, J, which for the 
case of proportionate coupling to the leads, r^(u;) = ar^(u;), can be expressed&iii 
in terms of the density of states, Po-(co'), as 

tl ^ J —CO 

where T„{lj) = r^(ci;)r^(ci;)/ [r^(cij) + r^(ci;)]. The density of states is given by 

p.(^) = --ImG;(^), (4) 

TT 

where G^{uj) is the Fourier transform of the retarded Green function, 

G:it) = -t9it){Mt),c:m). (s) 
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2. The Slave-Boson Hamiltonian 



Diagrammatic techniques are reliable when the expansion parameter is a small 
quantity. For an Anderson impurity with U — oo, it is natural to perturb in the hop- 
ping strength. However, the standard diagrammatic approach also requires that the 
unperturbed Hamiltonian be noninteracting, i.e., quadratic in the second-quantized 
operators. In the limit of infinite U, the bare Hamiltonian can be made quadratic 
by transforming the Hamiltonian (|T|) into a new Hamiltonian, expressed in terms 
of new local operators.!!!! These operators create the three possible states of the 
site: a boson operator b+, which creates an empty site, and two fermion operators, 
f+, which create the singly occupied states. The ordinary electron operators on the 
site, which transform the empty site into a singly occupied site or vice versa, are 
decomposed into a boson operator and a fermion operator, 

c.(t) =b+(t)f.(t) 

c:(t) = f;(t)b(t). (6) 

The "slave boson" in (H) acts as a bookkeeping device which prevents double occu- 
pancy of the site: when an electron creation operator acts on an occupied site, the 
boson part acting on the vacuum annihilates the state, c+ fi" \Q) = f+ bf^ \Q) = 0. 
(In these expressions \ fl) is the vacuum state.) 

In the "slave-boson" representation, the Hamiltonian for the infinite-?/ Anderson 
model becomes 

H= E e^^clc, + Y: e.f;f. + E (Vkacl^^i. + h.c). (7) 

The first two terms form the unperturbed, quadratic Hamiltonian and the last term, 
which represents hopping between site and leads, can be handled as a perturbation. 
The fermions and boson are treated as ordinary particles in the perturbation expan- 
sion. For example, the lowest order diagrams are shown for the boson and fermion 
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propagators in Fig. 0. While summation of a few low-order diagrams is possible,t§E3 
techniques are also available to sum whole classes of diagrams. In the end, whatever 
approach is taken, properties of the physical electrons can be constructed from the 
results for the boson and fermions. 

There is, however, an added constraint, as the site can only be in one of the 
states b+ \Q) and Accordingly, in all physical states, the number of bosons 

plus the number of fermions, 

Q = b+b + ^Cf., (8) 

a 

must be equal to unity. We will show below, when we describe the Keldysh dia- 
grammatic approach, how this constraint is dealt with. 



3. The Keldysh Formalism 

Previous diagrammatic calculations using the "slave-boson" represent at ion0 
have addressed equilibrium properties of the Anderson model. Since our focus is on 
nonequilibrium properties a different approach is required. Specifically, we employ 
the Keldysh§li rather than the Matsubara@ formalism. The main complication 
with nonequilibrium is that the basis of the equilibrium diagrammatic approach, 
the fact that the state of the system at t = +oo is identical to the state of the 
system att = — oo, uptoa phase (Gell-Mann and Low theoremil), is no longer 
valid. Since in a nonequilibrium system real dissipation can occur, the state of the 
system is in general not known at t = +oo and one must relate all quantities to 
the state of the system at t = — oo. In practice, this means that instead of having 
integrals from t = — oo to t = +oo as in the usual zero-temperature formulation,^ 
all integrals have to be carried out along a path that starts and ends at t = —oo 
(Fig. ^). Consequently, a Green function will depend not only on the times at which 
the operators act, but also on the corresponding branch of the contour. Thus the 
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Green functions carry additional indices, and the usual perturbation expansion, or 
the Dyson equation, takes a matrix form. In all, there are three independent types 
of two-particle Green functions in nonequilibrium. It is convenient to choose, in 
addition to the retarded Green function (|^), the two Green functions: 

G<{t) =z(cj(0)c.(t)) 

G>(t) = -z(c.(t)c+(0)), (9) 

as they carry information on the occupation of the site. 

For the problem at hand, the starting point at t = —oo is the Anderson impurity 
and the leads unconnected and separately at equilibrium, possibly with different 
chemical potentials. Formally, the hopping is turned on slowly, and nonequilibrium 
properties are evaluated long after the hopping is fully established, when a steady 
state has been achieved, but before current flow has changed the chemical potentials 
deep in the leads. 

Before applying the Keldysh formalism to the "slave-boson" Hamiltonian (0), 
one has to overcome the difficulty associated with the constraint that the physical 
states are restricted to the Q = 1 ensemble. Specifically, it is not convenient to 
perform diagrammatic calculations (Keldysh or Matsubara) in a restricted ensemble 
since Dyson's equation does not apply. Instead, we introduce a complex chemical 
potential,^ calculate diagrams in an unrestricted ensemble treating the hopping as 
a perturbation, and finally use the complex chemical potential to project to the 
Q = I subspace. In practice, this projection corresponds to keeping only an easily 
identified subset of diagrams. 

To formulate the Keldysh diagrammatic theory in terms of a complex chemical 
potential, it is convenient to start with a formal expression for expectation values 
in nonequilibrium. In the Q = 1 ensemble, the expectation value of an operator O, 
can be written as 

(0)q=i = ^Tr{e-^(^«-^-^--^-^-UQ,iTc[5c(-oo,-oo)0]}, (10) 
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where Tc orders operators along the Keldysh contour (Fig. ^ and the partition 
function is given by 

= Tr{e-'3(^«-^-^-~^«^-)(5Q,„Tc[5c(-oo,-oo)]}, (11) 



with Q = 1. In (10), the system evolves under the action of the Hamiltonian so that 

Sci-oo, -oo) = exp[-i <f dt'Hit')]. (12) 

J c 

Importantly, the operator O may include parts acting at different times, e.g., O = 
2C+(0) Ccr{t) would give the nonequilibrium expectation value of the "lesser" Green 
function, G'^{t). Since the Hamiltonian commutes with the sum of bosons and 
fermions, Q, the projection to the Q = 1 ensemble is accomplished once and for all 
by the factor 6q^i in (|10|). It is not necessary to include a chemical potential for 
the impurity since local expectation values in the coupled system are independent 
of the initial state of the impurity. 

To transform to an ensemble where Q is unconstrained,^ one rewrites the Kro- 
necker delta as an integral over a complex chemical potential, iX, 



5q, = ^ r^^dXe-^f^'^Q-'K (13) 



Dividing both numerator and denominator of ( pUj ) by Zq=q, gives^ 



Zq=\ 2tT J-n/l3 



|^(0)«, (14) 



where 



(0)a = ^Tr !e-^iHo-^^^N^-^.nNn+^XQ) Tc[Sc{-oo, -oo)0]\. (15) 

In (p^, {0)l]^ is the coefficient of the term of order exp{—ij3X) in {0)ix. The 
important point is that {0)i\ in (|T3|) is in the standard form for diagrammatic 
perturbation theory since the trace in (^) and in the partition function. 



= Tr {e-/5(^o-'^^^^-'^«^«+^^Q) Tc[^c(-oo, -oo)]}, (16) 

are taken over all states without restriction to Q = 1. 

According to ([l^) , the nonequilibrium expectation value of an operator, O, in 
the Q = 1 ensemble has two contributions: a normalization factor Zq=q/Zq=i, 
and the coefficient of exp(— z/5A) for the same operator in the iX ensemble. The 
normalization can be obtained from the identity ((5)q=i = 1, which implies 

|^ = (b^b)Sp + E(f;f.)Sp. (17) 

The expectation value {0)\]^ , as well as the expectation values appearing on the 
right-hand side of Eq. (^) can be obtained diagrammatically. 

B. Non-crossing Approximation 

To obtain a well-behaved density of states from the nonequilibrium perturbation 
theory, one needs some way of summing diagrams to all orders in the hopping. 
In finite-order perturbation theory there are divergences associated with the bare- 
levels, e„, and, at T = 0, logarithmic divergences near the chemical potentials due 
to the Kondo effect .0 To control these divergences, we employ the non-crossing 
approximation, which has been used successfully to treat the infinite-?/ Anderson 
model in equilibrium.i2l As can be seen from Fig. at lowest order in perturbation 
theory the boson self-energy involves the fermion propagator while the fermion self- 
energy involves the boson propagator. By using the two relations self-consistently 
- the non-crossing approximation, see Fig. | - one obtains a set of coupled integral 
equations, which can be solved numerically. Solving these self-consistent equations 
corresponds to summing a subset of diagrams to all orders in the hopping matrix 
element V. It can be shownH that all diagrams of leading order in 1/A^, where N is 
the number of spin degrees of freedom, are included in this subset. Therefore, the 
non-crossing approximation is expected to be a quantitative approach in the limit of 
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large A^. Importantly, comparison to exact result^ has shown that the non-crossing 
approximation is also quantitative for the case = 2, of interest for quantum 
dots. Here, we generalize the non-crossing approximation to nonequilibrium. The 
equations will involve not only the retarded Green function, but also the lesser and 
greater ones, leading to slightly more complicated equations than at equilibrium. 

Since our goal is to calculate the nonequilibrium current, we will calculate first 
the density of states for the Anderson impurity, Pcr{^^), and obtain the current from 
Eq. d^). To find the density of states for U oo, we need the retarded Green 
function (|^) in the ensemble with complex chemical potential (0), 
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Within the non-crossing approximation, the retarded Green function is expressed in 
terms of the full propagators for boson and fermion as 

-ze{t) [D>{-t)G^^{t)-D<{-t)G}^{t)], (19) 

where 

D>{t)^ -^(b(t)b+(0))!? 
D<{t)^ -^(bnO)b(t)),« 

G%{t)^ ^(f;(0)f.(t))!p. (20) 

Equation (|I9|) is straightforward to obtain by decomposing the electron operators 
into boson and fermion operators (^ and then factorizing the boson and fermion 
parts. The latter step corresponds to a neglect of vertex corrections. Since each 
term in ( [TP| ) contains exactly one "lesser" operator, with a boson or fermion lowering 
operator acting directly to the right, the overall result is 0(exp(— i/^A)) as required. 
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Because the Hamiltonian is time independent, it is simplest to evaluate the boson 
and fermion Green functions in the frequency representation. The physical density 
of states is then given by 

pAuj) = ^ ^ J^J^' [D'i^') + u') - D<{u') G)Au: + • (21) 

The non-crossing approximation is represented diagrammatically in Fig. ^: the 
boson and fermion propagators are each assigned a single self-energy bubble (al- 
beit determined self-consistently) and the self-energies are iterated to all orders via 
Dyson's equation. Standard manipulation of the nonequilibrium Dyson equations 
then leads toil 

D^{uo)= D'-{uo)Y{>{uj)D''{uo) 

G%{uj)=G}^iuj)J:l{uj)G%{uj), (22) 
where, the self-energies are given by 

U>{uj) = -- du'Y: \V,^\' gU^' - CO) G>^{uj') 



^A^)= ^/ duj'j:\Vka\'9U^-io')D>iJ). (23) 



keL,R 



In ([23|), the small g's are the Green functions of electrons in the leads not coupled 
to the site,il 



9kA^) = - 2vrz [1 - /L{fl)(^)] S{uj- Eka) 

9kA^) = 27ri /i(fl)(^) 6{uj- Eka). (24) 

Several other relations are required to close the equations for the non-crossing 
approximation. The retarded Green functions for the boson and fermions in (^2]) 
are given byil 

1 



w - Wiuj) 



^ w - - 
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where the retarded self-energies are Hilbert transforms of the greater self-energies 



27r J-oo u — uj' + irj 



Y'(uj) = — du' ' . (26) 

The advanced Green functions D'^ and G°j^ in (p2D are complex conjugates of the 
retarded Green functions D'' and G^^. Eqs. (pB]) follow because, by definition, all 
retarded Green functions and self-energies can be written as a difference of greater 
and lesser functions, G""(t) = 6{t)[G^{t) — G^{t)]. In the i\ ensemble, the lesser 
functions for the boson and fermions are 0[exp(— i/5A)] and must be dropped from 
the retarded functions which are 0(1). One therefore has the useful relations: 

G}^iu) = 2tlmG}^iu;), (27) 
for the boson and fermion Green functions, and 

n>(w)= 2ilmn'^(a;) 

S^^(a7) = 2zImS^^^(^), (28) 

for the self-energies. 

The closed set of equations for the non-crossing approximation can be solved 
iteratively. In practice, we start with an initial guess for the greater boson Green 
function, D^{lj), calculate Sj^(ci;) for each spin by combining (^) and (^^, and 
use ( P3| ) to get G^j^{uj). The values for the greater fermion Green functions, from 
G^^{uj) = 2ilmG^jr^{uj), can then be used in a parallel way to obtain an improved 
D^{uj). This procedure is iterated to convergence. A similar procedure is then 
followed for the lesser Green functions. Following an initial guess for D^{uj), the 
fermion self-energies, Sj^(a;) are obtained from (^), and Gj^{u!) is determined from 
(p^). The steps are repeated for D'^{uj), and the process iterated to convergence. 
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Finally, the physical density of states, Paioj) is evaluated by the convolution of the 
boson and fermion Green functions in (pT]). 



C. Numerical Methods 

In this section, we describe in greater detail the numerical procedures we have 
used to solve the nonequilibrium, self-consistent equations for the non-crossing ap- 
proximation. Following the equilibrium workj^l'll we take the energy dependence of 
the coupling between the site and the leads (H) to be Lorentzian 

The finite width, W , reflects the finite bandwidth in the leads and is necessary to 
prevent ultraviolet divergence of the results.^ In principle, the bands in the leads 
can be centered at different energies, but the validity of Eq. (|^) for the current 
requires V^i^uj) = aT^^u), so we take 6^ = 6^ = 69 = 0, throughout. [An expression 
( |Al| ) for the current in the absence of this condition is given in the appendix.] The 
choice of a Lorentzian form allows a simplification of the self-consistent equations.ii 
In general, to iterate the non-crossing approximation equations, the retarded self- 
energies for the boson and fermions (|26|) must be evaluated by double integrals over 
the greater Green function of the other species. For the the Lorentzian coupling, 
however, one of these integrals can be performed analytically. First, combining Eqs. 
(^) and (|26D , the boson retarded self-energy can be written as a single integral, 

/oo 
du;'H-^^-\uj'-u;)G%{u;'), (30) 

with the kernels 

+ ^ReUC-- : ^-^-^ ) - ^ ( ^ + " '^'r " ^^^^^^ ) j (31) 



■ I 



W 



27r 
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In (|3T|), [3 is the inverse temperature and '^{z) is the Digamma function.LZl Second, 
the fermion retarded self-energies can be written as single integrals 



/oo 
du:'K:;^-\u-u')D>{u:'), 



(32) 



with the kernels 



K. 



MR), 



(27r)2 



|vr[l (^)] 



2ti 



tie 



L(R) 



271 



(33) 



UJ 



w 



TT 

2 



)} 



2 27r 

Since the greater Green functions and self-energies are just the imaginary parts of 
the corresponding retarded functions (|27| - |28D , the above equations, together with the 
relation (^) between the retarded Green functions and self-energies, form a closed 
set. In practice, we make an initial guess for the greater boson Green function and 
then iterate the equations to convergence. Typically the results converge within five 
iterations. We have checked the accuracy of the results by comparing to the sum 
rules on the boson and fermion Green functions 



duo 



1 



7r 



duj 



--ImG},(^) 



(34) 



These relations are always satisfied to better than 0.5% by the converged numerical 
solutions. 

A separate iterative loop is required to evaluate the lesser Green functions and 
self-energies. Eqs. (|23D for the boson and fermion lesser self-energies can be rewritten 

as 



1 



(35) 



and 
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1 

-Y. du'Vr\uj-u')f,,,,{uj-J)D<{u'). (36) 

^TT J „ J—oo 



Together with Eqs. (^2]), these form a closed set of equations for the lesser Green 
functions and self-energies. Again, following an initial guess for the boson lesser 
Green function, these equations are iterated to convergence. Since the lesser Green 
functions have an arbitrary overall normalization, to check the convergence it is 
necessary to monitor a normalized quantity. We choose, for simplicity, to monitor 
the occupation of each spin state 

(n.) = |^[-f rdujG<^{uj)], (37) 

where the normalization is provided by the ratio of partition functions, which from 
Eq. (|17D is given by 

Zq=i i 



duo 



D<iuj)-j:GU^) . (3^ 



Zq=o 27r 

Typically, within five iterations the occupations converge to better than 0.01%. 
However, one cannot expect the accuracy of the results to be better than the accu- 
racy of 0.5% found for the retarded Green functions, from which the lesser functions 
are constructed via (P^). The final accuracy is verified by the sum-rule for infinite-?/ 
relating the total density of states of one spin state to the occupancy of all the other 
spin states 

/oo 
dujp^{u;) = l- (39) 

This relation is always satisfied to within 0.5%. 

The procedure has also been checked by comparing to the equilibrium results,!^ 
and excellent agreement is found. This is an independent check since much of our 
numerical procedure differs from that used in equilibrium. Importantly, because 
there are multiple sharp features in the nonequilibrium density of states (discussed 
in the following section), we have used a self-adjusting mesh for the numerical in- 
tegrations rather than the logarithmic mesh used in the equilibrium case.i3ii We 
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have also found that evaluating the Green functions on the range [— lOW, lOW] is 
sufficient for numerical accuracy. The wider range used previouslyH'll is unneces- 
sary because the kernels ( pT]) and (^) contain all the effects of the long band tails 
in the leads. 



D. Results of the Non-crossing Approximation 

In this section, we present numerical results of the non-crossing approximation 
for an Anderson impurity in and out of equilibrium. 



1. Linear Response Conductance 

The equilibrium properties, calculated by the non-crossing approximation, can 
be used to predict quantitatively the lineshape of the linear-response conductance 
peak. To our knowledge, this is the first detailed prediction of the conductance peak 
evolution in the Kondo regime. From Eq. (|^) the linear-response conductance is 
given by 

^ = 27r-r E / duj [-f'{w)] T^{uj) p.(^), (40) 

where Po-(cu) is calculated at equilibrium, and j'iw) is the derivative of the equi- 
librium Fermi function. In Fig. |] we plot the equilibrium density of states (|21|) of 
an Anderson impurity with two degenerate spin states, for one value of chemical 
potential (/i — eo = 2r). There is a sharp Kondo peak at the chemical potential. Its 
amplitude increases with decreasing temperature down to the Kondo temperature,^ 
~ iy(r/27r(/i - eo))^/^exp[-7r(/i - eo)/r], where it saturates. In Fig. |(a), the 
resulting linear-response conductance obtained from ( ^OD is plotted as the chemical 
potential is swept through the bare-level energy at three different temperatures. Sev- 
eral features are noteworthy. First, as the temperature is initially lowered, the width 
of the conductance peak decreases proportional to /cgT, because the peak hneshape 
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is determined by the derivative of the Fermi function ( pOD . As the temperature is 
lowered below F, the peak width is expected to be dominated by F and saturate. 
Here, however, for kgT ~ 0.075F, the conductance peak begins to broaden again. 
This broadening is entirely due to the appearance of the Kondo peak in the density 
of states, and is therefore a signature of the Kondo effect. Second, as tempera- 
ture is lowered, the peak amplitude increases and finally saturates. The saturated 
open channel conductance, a = 2e^//i, is only achieved for a dot symmetrically cou- 
pling to its leads; otherwise, the conductance is reduced by the asymmetry factor 
4Fq^FJ^/ (Fq^ + F(^^). Third, the peak maximum shifts to higher chemical potentialS 
and the tails become power law (roughly Lorentzian) rather than exponential as at 
higher temperatures. It is interesting to note that the conductance peak remains 
nearly symmetric with decreasing temperature despite the very asymmetric behavior 
of the density of states, which has a Kondo peak only for fi > cq. For comparison, 
the total occupancy of the site is plotted as a function of chemical potential in Fig. 
^(b). Unlike the conductance, the total occupancy is not sensitive to the behavior 
near the Fermi surface, and therefore does not show any obvious signature of the 
Kondo effect. 

The temperature dependences of the main features of the conductance peak are 
plotted in Fig. 0. Over a broad range of temperatures the peak width, amplitude, 
and position increase roughly logarithmically with temperature. This reflects the 
logarithmic scaling of interactions which is the well known signature of the Kondo 
effect in perturbation theory.!! For kgT ~ 0.005F, the peak amplitude saturates 
while the peak position and width continue to increase. Note that the non-crossing 
approximation is known to overestimate the Kondo peak amplitude somewhat for 
chemical potentials within a few F of the bare-level energies.ii The true magnitude 
of the conductance peak for a symmetric structure is therefore not expected to 
approach the maximum value, 2e^//i, until temperatures below those shown in Fig. 
p. For example, taking the saturated low temperature value for (rio-) from Fig. 



P(b), which is known to be rehable,L3 and using Langreth's exact relation for zero 



one obtains a conductance a = l.GSe^/h at the peak of the lowest temperature 
curve in Fig. ^(a). The observed peak value of a = l.dOe^/h is therefore 15% higher 
than the expected zero-temperature value, and must be an overestimate. However, 
the tendency of the non-crossing approximation to overestimate the conductance 
vanishes as 2{n„) — > 1. 



There are qualitatively new features in the nonequilibrium density of states com- 
pared to equihbrium. In Fig. |], the density of states (^) of an Anderson impurity 
with two degenerate spin states is plotted both for equilibrium and for nonequilib- 
rium, where the two leads have different chemical potentials. There are striking 
differences between equilibrium (solid curve) and nonequilibrium (dashed curve). In 
equilibrium there is a single Kondo peak at the chemical potential. Out of equilib- 
rium, the Kondo peak splits into two smaller peaks one at each chemical potential. 
With decreasing temperature, the amplitudes of these peaks do not increase to the 
unitarity limit, but saturate at a much lower value. This saturation occurs at a 
temperature above 71-, and results from dissipative processes in which an electron 
is transferred from the lead with higher chemical potential to the lead with lower 
chemical potential. The non-crossing approximation includes these processes since 
it has contributions from all orders in the hopping, but separating out the relevant 
diagrams is not straighforward. Instead, in the next section, we present an analytical 
formula for the dissipative lifetime obtained via perturbation theory. 

A clear signature of the Kondo effect is expected in the nonlinear current. For 
chemical potentials above the bare-level energy, the linear-response conductance 





(41) 



2. Nonequilibrium 
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is dominated by the narrow Kondo peak in the density of states. In nonhnear 
response, at low temperatures, the current, J, is determined by an integral of the 
density of states between the two chemical potentials (^. Therefore, as soon as the 
chemical potential difference exceeds the width of the Kondo peak, the differential 
conductance will fall off dramatically. Moreover, the Kondo peak will split and 
the split peaks decrease in amplitude with increasing chemical potential (Fig. |^). 
The net effect is a sharp maximum peak in the differential conductance around 
zero bias, for /i > eo.&S In Fig. ||, we have plotted the differential conductance 
as a function of chemical potential difference, or equivalently voltage bias at two 
temperatures. The expected peak is clearly resolved at a temperature kgT ~ O.OSF. 
This is substantially higher than the temperature, k^T ^ 0.025F, at which the 
linear-response conductance peak has broadened unambiguously (+10%) over the 
minimum width. The peak in the nonlinear differential conductance is therefore 
likely to be the first signal of the Kondo effect in transport through a quantum dot. 

It is worth noting that the observability of the Kondo peak in the differential 
conductance depends only on the ratio kgT/T, not on the Kondo temperature, 2]^. 
To demonstrate this, the differential conductance is plotted at different tempera- 
tures in Figs. ^ and for level depths differing by F, and, consequently, Kondo 
temperatures differing by a factor of exp(7r| Aeo|/F) ^ 23. If the presence of a zero 
bias-peak depended on the Kondo temperature one would expect the peaks to wash 
out at temperatures differing by a factor of 23. Instead, both the zero-bias peaks 
are clearly visible at k^T = 0.05F (Fig. g), but by k^T = O.IF (Fig. |10D both 
peaks have washed out. This behavior of the differential conductance reflects the 
temperature dependence of the Kondo peak in the density of states. While the final 
saturation of the density of states peak occurs at temperatures below the Kondo 
temperature, ~ exp(— 7r(yU — eo)/F), the temperature at which the Kondo peak 
first appears depends only on the coupling strength F. Since the peak first appears 
below kgT ~ O.IF, the differential conductance develops a zero-bias peak just below 
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this temperature. 



III. DISCUSSION 
A. Theoretical Results 

The most important result of this paper is a better quahtative and quantitative 
understanding of the low-temperature nonequilibrium properties of an Anderson 
impurity. The nonequilibrium characteristics, in particular the transport properties, 
follow from the form of the nonequilibrium density of states, Po-(co'). In this section, 
we discuss the main features of the density of states using both the results of the 
non-crossing approximation and other methods.&iil 

The most obvious features in the low-temperature equilibrium density of states 
for an Anderson impurity with /i > eo are the sharp peak at the chemical potential 
and the low, broad peak around the bare-level energy (Fig. To understand these 
features it is useful to recall how, at equilibrium, the density of states depends on the 
eigenstates of the system. At T = 0, the density of states, Pcr(t^) = —(1/^) I^^aX^); 
involves transitions from the A^-particle ground state to all + 1 or — 1 particle 
states. Since the correlated ground state of an Anderson impurity has a finite 
amplitude to have an empty site, the density of states includes a narrow peak due to 
transitions from the A^-particle ground state to the ground state with one more or one 
less electron. By definition the ground-state energies differ by the chemical potential, 
so this Kondo peak in the density of states occurs at the chemical potential. The 
weight of the Kondo peak is small, however, since the probability that the site 
is unoccupied in the ground state is much less than one. The remaining weight, 
associated with transitions to excited states, forms the low, broad peak around the 
bare-level energy, eo. For finite interaction energy, U, there is an additional broad 
feature in the density of states near eo -l- U; this feature does not appear in Fig. |^ 



21 



because of the limit f/ — > oo. 

Out of equilibrium, there is no true ground state of the system, but quantum 
fluctuations still produce a finite probability of an empty site. As at equilibrium, 
these fluctuations involve electrons hopping between the site and states in the leads 
near each chemical potential.i'i The N ^ N + 1, and ^ — 1 transitions which 
determine the nonequilibrium density of states therefore include some excitations 
which change the system only by adding an electron or hole near one of the chem- 
ical potentials. These low-energy transitions produce the Kondo peaks near each 
chemical potential in the nonequilibrium density of states. 

Unlike equilibrium, the configurations of the system out of equilibrium are not 
true eigenstates, but have a finite lifetime, r. The energies of transitions are there- 
fore broadened by h/r, and all features in the density of states are broadened an 
equivalent amount. This is the origin of the suppression of the Kondo peaks out of 
equilibrium (Fig. |^). The finite lifetime is due to real processes in which an electron 
is transferred from the higher chemical potential lead to the lower chemical potential 
one. An estimate of this lifetime can be obtained from straightforward perturbation 
theory in the coupling strength. One assumes that the site is initially occupied by an 
electron of spin a and calculates the decay rate using the Golden Rule.0 The only 
complication is that, at T = 0, the lowest order, energy-conserving process involves 
two separate tunneling events (the site electron hops out, and another electron hops 
in), and therefore occurs at O(V^). Allowing for two, possibly nondegenerate spin 

states, we find 
1 1 ^ „ 1 

^ A,B=L,R 



^ 1 1 roo 

r = fi Er;:'(^<^)[i-/AM] + i-r E / 



1 1 

+ 



(e^ - e + iT]y (e^ - e - iijY 



r^(e) rf,(e - + e^O [1 - /^(e)] /^(e - + t,, 

(42) 

For a deep level, cq < HdI^ri so at zero temperature and for constant F, Eq. (H2p 



reduces to 
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- = ^ E r^r^.'dif^B-f^A + e.-e.,) ^^ + \ (43) 

which exphcitly shows that the hfetime is only nonzero for finite bias, or finite 
level splitting. The results of the non-crossing approximation are consistent with a 
broadening of Kondo peaks by the inverse of the nonequilibrium lifetime h/r^. In 
effect, h/r^ is an new cutoff energy for the logarithmic scaling of interactions in the 
Kondo problem.ii 



B. Relation to Experiment 

Since a quantum dot weakly coupled to its leads is an Anderson impurity, the 
results of the previous sections have practical significance. Specifically, at sufficiently 
low temperatures, transport through a quantum dot will be dominated by the Kondo 
effect. We discuss the practical requirements for the Kondo effect to be observed in 
quantum dots, and suggest possible experiments. 

There are two general classes of transport experiments to study the Kondo ef- 
fect in quantum dots: linear response and nonlinear response. While detection of 
the Kondo peak in the density of states is possible in the linear-response conduc- 
tance, the nonlinear conductance offers a clearer signature and one that persists to 
higher temperatures. Figs. ^ and ^ indicate the appearance of the Kondo effect 
in linear response. The sweep of chemical potential indicated in the figures can be 
accomplished by sweeping the voltage of a separate gate which couples capacitively 
to the dot.lii In fact, conductance peaks with Lorentzian tails have already been 
observed in transport through a quantum dot by this method.ill The long tails of 
the peaks imply that coherent transport of electrons is taking place. However, so 
far no broadening of the conductance peaks at zero magnetic field is observed down 
to T ~ 50mK. This is consistent with the prediction of the non-crossing approxima- 
tion that noticeable broadening (+10%) occurs by kgT ~ 0.025r, since the largest 
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resonance width in the experiment is F ~ 40/ieV for which 10% broadening is not 
reached until T ~ lOmK. 

The appearance of the Kondo effect in nonhnear response is shown in Figs, 
^, and |10|. The sharp drop of the differential conductance around zero applied bias 
reflects the sharpness of the Kondo peak in the density of states.&S Furthermore, the 
peak in the differential conductance persists to kgT ^ 0.05F and therefore should be 
observable in existing quantum dots up to T ^ 20mK. The magnitude of this zero- 
bias peak is optimized by performing the differential conductance measurement at 
the half-maximum point of the linear-response conductance peak (as we have done 
in Fig. H). 

An additional, striking signature of the Kondo effect in nonlinear response is 
the evolution of the peak in the differential conductance with magnetic field. From 
perturbation theory,@ and from an equations-of-motion approach,&i it can be shown 
that a finite magnetic field shifts the Kondo peaks in the nonequilibrium densities 
of states by the Zeeman energy, and consequently splits the peak in the differential 
conductance by twice the Zeeman splitting of the levels. 

Unfortunately, the behavior of the differential conductance in a finite magnetic 
field is beyond the scope of the non-crossing approximation. Specifically, when the 
level degeneracy is broken, the non-crossing approximation produces, in addition 
to the peaks found by the other methods, spurious peaks in the density of states. 
These peaks are due to a false interaction of each level with itself, brought on by the 
neglect of vertex corrections. While self-interaction effects are unimportant in the 
large- iV limit, for finite the corrections can be significant. An extreme example is 
the non- interacting case, A^ = 1, where the non-crossing approximation incorrectly 
predicts a Kondo peak in the density of states. Because of this false self-interaction, 
the non-crossing approximation produces additional Kondo peaks at the chemical 
potentials for non-degenerate levels, and therefore is unreliable for transport prop- 
erties in a magnetic field. Interestingly, the non-crossing approximation continues 
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to produce reliable results in a magnetic field for thermodynamic quantities {e.g. 
magnet iz at ionil) which depend on the entire density of states and not just on the 
behavior near the Fermi surface. 

IV. CONCLUSION 

In summary, we have analyzed the low-temperature, nonequilibrium properties 
of an Anderson impurity in the limit of infinite on-site interaction. The model 
corresponds to a quantum dot, weakly coupled to two leads with different chemical 
potentials. The Kondo effect, which dominates transport through the impurity, is 
modified by two new energies present in nonequilibrium: the chemical potential 
difference A/i, and the inverse of the dissipative lifetime h/r^ (|42|). These energies 
are apparent in the nonequilibrium density of states, which we obtain via the non- 
crossing approximation. The chemical potential difference appears in the density 
of states via the splitting of the Kondo peak into two peaks, one at each chemical 
potential. The amplitudes of these peaks are suppressed by dissipative processes in 
which an electron is transferred from the higher chemical potential lead to the lower 
chemical potential one. 

Experimentally, we predict that the Kondo effect can be observed in transport 
through a quantum dot by either linear or nonlinear measurements. The emer- 
gence of the Kondo peak in the density of states at low temperatures will cause the 
linear-response conductance peaks vs. gate voltage to broaden, shift, and increase 
in amplitude roughly logarithmically with decreasing temperature. For a symmetric 
structure the conductance amplitude will saturate at 2e^//i, the conductance of an 
open channel.0lli The clearest feature of the Kondo effect in linear response, how- 
ever, is the broadening of the conductance peak, which is predicted to reach +10% 
below kgT ~ 0.025r, where F is the total coupling strength to the leads. In non- 
linear response, the Kondo peak will produce a peak in the differential conductance 
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around zero bias.Q Since this nonlinear peak remains clearly defined for temperatures 
up to kgT ~ O.OSr, we believe it will be the most accessible signature of the Kondo 
effect in quantum dots. 

We thank Dan Cox, Roger Lake, and Patrick A. Lee for valuable discussions. 
Work at U.C.S.B. was supported by NSF Grant No. NSF-DMR90-01502 and by 
the NSF Science and Technology Center for Quantized Electronic Structures, Grant 
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APPENDIX A: 



In this appendix, we show that the non-crossing approximation is current con- 
serving. Specifically, the current through the Anderson impurity can be expressed 
either as a current flowing from the left lead into the site or as a current flowing 
from the site into the right lead. Within the non-crossing approximation, these two 
expressions are equivalent. 

The full expression for the current through the left (right) tunneling barrier, with 
no restriction on the relative couplings to the leads,0 is 

< 

where G^{uj) are Fourier transforms of the physical-electron Green functions deflned 
in (^). In the non-crossing approximation, these Green functions factorize into boson 
and fermion parts, 

atico) f 1^ rdu'D<{uo') G%{co + CO'). (A2) 

Since this factorization is an approximation, one can ask whether the two expressions 



for the current (|A1|) remain identical. 

To show that current is conserved in the non-crossing approximation, we examine 
the difference, — J^, between the currents flowing through the two tunneling 
barriers, 

Ann ZjQ = \ o- ^^J-OO J-OO ^^g-J 

X I [l - j,U^)\ D>{uo') G)X^ + CO') + f,U^) D<{uj') G%{co + co')^ 

This expression is simplifled by recognizing that the cj-integration produces factors 

of the boson self-energies, 
1 

n>(^')= ^EE/ rf^rrM/.(.)MG^^.(^ + ^') 



a L,R ' 



1 ,00 , , (A4) 

n<(o.')= -^EE/ rf^r^H i-/.(.)M + 
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The difference can therefore be written as 



J L Jr 



1 Zc 



duj' 



D>(uj')U<(u')-D<(uj') U>(uj') 



(A5) 



2n Zq=i 

The integrand vanishes because of the relation between the boson Green functions 
and self-energies (P^l), 



D>M = W{uj) D\uj) 



(A6) 



The non-crossing approximation therefore explicitly conserves current, 



J T. — J R 



(A7) 
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FIGURES 

FIG. 1. Schematic band diagram of a quantum dot coupled via tunneling barriers to two leads 
with different chemical potentials. At zero magnetic field, the energy level eo on the quantum dot 
will be spin degenerate, and a large Coulomb interaction energy, U, will prevent double occupancy. 

FIG. 2. Diagrammatic expansion for (a) the slave boson and (b) the fermion propagators. 
The coupling between site and leads is treated as the perturbation, so each vertex corresponds to 
a tunneling event. 

FIG. 3. Real-time contour for nonequilibrium Green functions in the Keldysh formalism. 

FIG. 4. Diagrammatic representation of the non-crossing approximation, (a) Dyson's equation 
for the boson propagator includes the fermion propagators in the self-energy, and (b) Dyson's 
equation for each fermion propagator includes the boson propagator in the self-energy. 

FIG. 5. Equilibrium and nonequilibrium density of states for an Anderson impurity 

symmetrically coupled to two leads of Lorentzian bandwidth 2W and chemical potentials /Xj, and 
flu- The impurity has two degenerate spin states at energy eo = 0, and an on-site interaction 
U ^ oc. With all energies in units of the total coupling to the leads, F, the band half- width is 
W = 100 and the temperature is T = 0.005. At equilibrium (solid curve) there is a single Kondo 
peak in the density of states at the chemical potential = = 2 (see inset) . Out of equilibrium 
(dashed curve), the peak splits into two suppressed peaks, one at each chemical potential, = 2.4 
and //fl = 2. 

FIG. 6. (a) Linear-response conductance a through an Anderson impurity for three different 
temperatures as a function of chemical potential. The impurity has two degenerate spin states at 
eo = 0. The conductance peak first narrows then broadens with decreasing temperature, (b) Total 

site occupancy, + n^, as a function of chemical potential for the same temperatures. 



29 



FIG. 7. (a) Temperature dependence of linear-response conductance peak position, (b) Tem- 
perature dependence of conductance peak amplitude, (c) Temperature dependence of conductance 
peak full-width at half maximum. In all three panels, the non-crossing approximation results are 
the data points and the solid curve is a guide to the eye. For comparison, the dashed curves are 
the exact results for non-interacting levels. 

FIG. 8. Differential conductance, edJ/dA/x, with /j.^ = 1.9, vs. applied bias, at two tempera- 
tures, ksT = 0.005 and fc^T = 0.05. The peak in the differential conductance at zero bias reflects 
the Kondo peak in the density of states. 

FIG. 9. Differential conductance, edJ/dA/j, vs. applied bias for /x^ = 1.9 (solid curve) and 
fiji = 2.9 (dashed curve) at kgT = 0.05. Zero-bias peaks due to the Kondo effect appear for both 
curves despite the very different Kondo temperatures. 

FIG. 10. Differential conductance, edJ/dA/i vs. applied bias for /x^j = 1.9 (solid curve) and 
fiji = 2.9 (dashed curve) at kgT = 0.1. Both zero-bias peaks in Fig. 9 have become shoulders. 
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